Effects of phonon interaction on the pairing in the high-Tc superconductors 
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We study the effects of phonon interaction on the superconducting pairing in the background of a d-wave 
gap, mediated by antiferromagnetic (AFM) spin fluctuations, using coupled BCS gap equations. We found that 
phonon interaction can induce a s-wave component to the d-wave gap in the (D+S) form with an interaction 
anisotropy and in the (D+iS) form without anisotropy, respectively. In either case, however, T c is not enhanced 
compared to the pure d-wave pairing without phonon interaction. On the other hand, anisotropic phonon interac- 
tion can dramatically enhance the d-wave pairing itself and therefore T c , together with the AFM spin fluctuation 
interaction. This (Dafm + D p h) type pairing exhibits strongly reduced isotope coefficient despite the large 
enhancement of T c by phonon interaction. Finally, we study the combined type of (Dafm + D p h +iS)) gap and 
calculate the penetration depth and specific heat to be compared with the experiments. 
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I. INTRODUCTION 

The problems of the high T c superconductors (HTSC) 
have been mostly focused on the electronic correlations and 
phonons were usually considered as a secondary player at 
best. However, recent Angle Resolved Photoemission Spec- 
troscopy (ARPES) experiments revived the interest of the 
possibly important role of phonons in the high-T c cuprates 
(HTC) 1 . In particular, the systematic measurements and 
analysis 2 of the kink structures in the quasiparticle disper- 
sion near the Fermi surface (FS) proved that: (1) the origin of 
the kink is electron-phonon coupling and its coupling is very 
strong, (2) the typical energy of the phonon(s) is ~ 40 - 70 
meV, and (3) the coupling matrix is quite anisotropic. There- 
fore, it is currently a pressing question what possible roles and 
effects, in particular, for the superconducting (SC) pairing, the 
phonon(s) can do in the HTC. 

Although the SC gap symmetry in the high-T c cuprates 
is well established as a d-wave by most experiments, there 
are also continuous experimental reports such as tunnelling 
conductance 3 , penetration depth measurements 4 , etc that pro- 
vide rather convincing evidences for a small magnitude of a s- 
wave component in addition to the dominant d-wave SC gap in 
some of HTC compounds. And it is well known that phonons 
mediate an attractive interaction between electrons and lead to 
an isotropic s-wave gap. Therefore, more sharpened question 
will be: is it possible to reconcile a d-wave gap and a s-wave 
gap together ? and in that case, what is the role of phonon(s) 
7 

While there is no consensus yet on the pairing mechanism 
for the d-wave gap in HTC, in this paper we assume that the 
AFM spin fluctuations is the mediating glue for d-wave gap 5 . 
The key ingredient of this mechanism is that the sign chang- 
ing character of the d-wave gap turns the all positive definite 
potential of the AFM spin fluctuations mediated interaction 
(we will call it "AFM interaction" for short from now), in mo- 
mentum space and in spin singlet channel, into an effectively 
attractive pairing interaction. From this point of view, it is 
clear that the phonon interaction and AFM interaction are an- 
tagonistic and difficult to work together for making a pairing 



because they tend to promote a different symmetry of the SC 
gap each other. 

An easy way to reconcile the d-wave and s-wave gap with 
two different pairing interactions would be to invoke two sep- 
arate bands for each gap 6 . The similar idea for the two s-wave 
gaps was very successful for MgB2 7 because there two well 
separated bands exist. However, in HTSC, numerous experi- 
ments, in particular ARPES measurements 8 , show that there 
is only one main band crossing FS. Then any idea of having 
s-wave and d-wave gaps in the single band, induced by two 
very different pairing interactions, seems too naive. How- 
ever, if the coupling matrix of phonon interaction possesses 
a strong anisotropy, this combined pairing problem is not a 
trivial one and needs systematic investigation with a traceable 
formulation. This problem has been alread y studi ed by many 
authors with different techniques and aim d 9 ! 10 ! 1 ^. In this pa- 
per, we took a simplest approach and studied the coupled gap 
equations for multiple gaps of a single band with two pairing 
interactions, i.e., AFM interaction and anisotropic phonon in- 
teraction, within the BCS framework. More justification for 
this approach will be discussed in the next section. By numer- 
ical solutions, we extensively investigated the necessary con- 
ditions of the strength and the degree of anisotropy of phonon 
interaction for various multigap solutions. 

We found that the (D+S) and (D+iS) type solutions are in- 
deed possible with a proper degree of anisotropy of phonon 
interaction for each case. However, T c is not enhanced at all 
in these cases even though an additional s-wave component is 
formed due to the phonon interaction. Other interesting pos- 
sibility is that the anisotropic phonon interaction can induce 
an additional d-wave component (D p h) to the AFM interac- 
tion induced d-wave component (Dafm)- As a result, T c is 
dramatically enhanced due to the phonon interaction. We then 
derived an analytic T c equation for this (Daf m + D p h) case 
and showed that phonon isotope effect can be strongly reduced 
due to the interplay between the AFM and phonon interactions 
despite the large increase of T c by phonon. This result can 
explain a long standing puzzle of the small phonon isotope 
effect in HTSC. In view of experiments^, the best anisotropic 
phonon can be the Bi ff buckling phonon mode which fulfills 



2 



all necessary qualifications for our model. Similar conclusion 
was obtained by other authors 10 using different approaches. 
We emphasize that the final effect of phonons in HTC com- 
pounds should be judged by an average effect of all important 
phonons together, which is beyond the scope of this paper. 

Finally, we also considered the combined type of (Dafm + 
D p h + iS) gap which is a very natural solution of the cou- 
pled gap equations and shows features such as (1) large en- 
hanced T c , and (2) appearance of a s-wave component at 
low temperatures. This type of gap will be a possible solu- 
tion for the tunnelling conductance experiment in YBCC 3 . 
Calculations of superfluidity density also demonstrated that 
(Dafm + D p h + iS) gap can explain the non-monotonic tem- 
perature dependence of the penetration depth measurements 
in YBCO and LSCO compound^. 



n. FORMALISM 

The interplay between phonons and electronic correlations 
in HTSC has been studied by numer ous aut hors using differ- 
ent theoretical approaches and model d™ 10 ! 11 !. In this paper, we 
took a simple minded approach with a specific question for 
this problem. Experiments tell us that, after all correlation 
effects taken place, the quasiparticles (q.p.) remain to form a 
fermi surface (FS) 8 . And quite detailed properties of the AFM 
spin fluctuations and symmetry allowed phonons are experi- 
mentally measured 1 ; theses measured bosonic fluctuations are 
final outcomes after all interplays between them and correla- 
tion effects. Now on the phenomenological basis, we consider 
a q.p. band and the AFM spin fluctuations as a dominant pair- 
ing interaction. In addition to that we add anisotropic phonon 
interactions and we study the pairing instability of the model 
for several possible types of multi-gap solution using a gen- 
eralized BCS gap equations. In this approach, still remaining 
interplay between spin fluctuations and phonons is ignored. 
The results of this paper should be taken with this caveat. The 
Hamiltonian is written as 



k <* kk'U 
kk'U 



and 



(1) 



where e(k) is the dispersion of the quasiparticles created by 
cjL as standard notation. Vaf m {k, k ) and V p h(k,k ) are 
the effective interactions, for the singlet superconducting pair- 
ing channel, originating from the AFM spin fluctuations and 
phonon(s), respectively. For traceable numerical calculations, 
we further simplify the above Hamiltonian as follows. The 
real two dimensional FS is simplified as a circular FS and the 
interactions are also modelled accordingly as follows. 



v / Aw a _ / ~Vp for \A(f>\ < 4>AN 

v ph {M) -| for m > ^ an 



(3) 



where A<f> =(/) — (/) and <Pafm — tt/2 representing the ex- 
changed momentum k k and the AFM ordering vector Q 
in the circular FS, respectively. Vm and Vp are chosen to be 
positive, so that Vafm is all repulsive and V p h is all attractive 
in momentum space. 

The property of the AFM interaction with a short range 
correlation is simulated with the inverse correlation length 
parameter <pQ (A = (na/y/2)(f>o). In this paper, we chose 
0o=l (A ~ 2a) for all numerical calculations, which is quite 
a short range AFM correlation. The degree of anisotropy of 
the phonon interaction is controlled by the anisotropy angle 
parameter 4>an which restricts the scattering angle between 
incoming and outgoing q.p momenta; for example, c/)an = ^ 
would allow a perfectly isotropic interaction. The phonon in- 
teraction V p h{4> — <f> ), however, does not restrict the incom- 
ing ((/>) and outgoing momenta (</> ). Now the reduced BCS 
Hamiltonian in the mean field theory can be written as 



H = ^40 c lt; < j c <t>'i<? + ^2 A AFM(<t ) )c<pa c -<Pa 

+ Y1 A * P h( ( t ) ) c <i>a c -<Pa 
<t>t 



(4) 



Vafm(A^) = V M 



(A0±. 



'AFM 



(2) 



where Aafm(4>) is the SC gap function induced by Vafm 
and A p h((p) is the one induced by V p h- The two gap functions 
a afm(4>) and A p ^(0) may or may not have the same sym- 
metry. After diagonalizing the above Hamiltonian we obtain 
two self-consistent equations as 

A A fm{4>) = Y Vafm(4>- (t> ) < c^|C_0 5T >, (5) 
A p/iO) = ^2v ph (4>- <j>) < c 0a c_ 0CT > . (6) 

Similar coupled gap equations were studied in previous 
studies 6 . The key difference is that our model has only one 
band. Then because the same band electrons should form the 
multigaps, there is severe competition between different gaps 
and more constraint to allow multigap solutions. In the next 
sections, we will consider several possible multigap solutions 
of the above coupled gap equations. 



A. (D+S) case 



In this case, we assume Aafm(4>) — AdCOs(2</>) and 
A p ft,(0) = A s . This leads to the two coupled gap equations 
as follows. 
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FIG. 1 : (Color online) (D+S) case. Calculated magnitudes of d-wave 
components (Ad/cjafm) an d s-wave components (As /uiafm) for 
different anisotropy angles 4>an = 7r/4, vr/2, and 7r. For all cases 
AT(0)Vm = 2.0, 7V(0)Vp = 0.45 and u ph /uj A FM=Q-5. 



J2Vvh(<f>-<i>')M4>'M<t>',u P ), (8) 



%(0 .WAFM,^) = N(0) 



•*» tegh(^) 



where £(</>) = ^/f 2 + A? 0), and A t (0) = (A d cos(2</.) + 
A s ) is the total gap function. iV(0) is the density of states 
(DOS) at FS and u>AFM,ph are the BCS energy cutoffs of each 
pairing interactions Vafm and V p h, respectively. Nonzero 
value of A s reduces the symmetry S4 of the pure d-wave gap 
to C2; namely, the nodal points shift away from diagonal di- 
rections and the sizes of the positive lobe and negative lobe of 
the total gap function A t (<j>) become different. We found that 
the coexistence of d- and s-wave gaps is possible only when 
V p h{4> — 4> ) is strongly anisotropic; the minimum anisotropy 
of the phonon interaction is 4>an < 7r/2 for our model inter- 
action of Eq.(3). 

Fig.l shows typical results of the (D+S) gap solution. 
The normalized d-wave gap Ad /u>afm and the s-wave gap 
A 5 /ujafm are plotted for different values of anisotropy angle 
4>an — 7r/4, 7r/2 and tt with fixed values of coupling con- 
stants N(0)V M = 2 and N(0)V P = 0.45. The effective 
dimensionless coupling constants Xafm,pH are smaller than 
theses values which is defined below as a projected average 
with a SC gap function. 



where r)((f>) = cos(2</>) for d-wave gap and r)(<p) — 1 for s- 
wave gap. For the chosen potential strengths N(0)Vm = 2 
and N(0)Vp = 0.45 in Fig.l, the effective coupling strengths 
are \ A fm,d = 0.332 and \ phtS = 0.450, 0.338, 0.198 for 
4>an — T, 7r/2, 7r/4, respectively. 

When 4>an = 1" (corresponding to a perfectly isotropic 
phonon interaction), solution is either a pure d-wave or a pure 
s-wave gap primarily depending on the strength of Xafm and 
\ p h\ the final dominant instability is determined not only by 
the interaction strengths but also with the cutoff energy scales 
Uafm and uj ph . In this case, \ p h,s (=0.450) is much bigger 
than Xafm,d (=0.332), therefore a pure s-wave gap becomes 
a solution. Only when <Pan < tt/2 and X p h,s is not dominant 
over Xafm,d, a s-wave gap component can coexist with a d- 
wave component. When these conditions satisfied, no matter 
how weak the phonon interaction is, a finite s-wave compo- 
nent coexists over the entire temperatures and shares the same 
transition temperature T c . Also this T c remains the same as 
the pure d-wave T c o without the phonon interaction, regard- 
less of the presence of As and its magnitude. This some- 
what unexpected result can be understood from the coupled 
gap equations Eq.(7-8). As far as the d-wave component is 
nonzero, that is the dominant pairing gap and T c is deter- 
mined solely by the d-wave gap equation Eq.(7) in the limit 
of Ap, As — > 0. In the opposite limit, namely, when X p h.s is 
dominant over Xafm.d, then even if <Pan < tt/2, the d-wave 
gap is entirely suppressed and T c is determined solely by the 
s-wave gap equation Eq.(8). 

In Fig. 2 we fixed the anisotropy angle 4>an{= tt/4) and 
(9) changed the coupling strength of phonon N(G)Vp (=0.7, 0.9, 
1.1). The variation of the magnitude of the s-wave gap As 
can be understood with the effective coupling strength of 
X p h,S (=0.31, 0.39, 0.48, respectively) compared to Xafm,d 
(=0.332) as in Fig.l. One peculiar feature occurs when 
N(0)V P = 0.9 (A p/l! s=0.39)(down triangles). The higher 
transition temperature T c ,high is understood as the coexis- 
tence cases of Fig.l; although X p h,s (= 0.39) > Xafm d (— 
0.332), the fact ujafm > Wph makes the d-wave pairing still 
slightly dominant at high temperatures. But at lower temper- 
ature one more transition occurs where the d-wave gap col- 
lapses to zero and the s-wave gap component makes a sudden 
increase. This abrupt change of gaps is a first order transition 
and indicates that there are two competing local minima in the 
free energy. With temperature the global minimum changes 
from one local minimum to the other local minimum. This 
peculiar behavior is due to the closeness between Xafm,d 
and X P h,s, and yet the distance of the energy scale between 
ojafm = 1 and up = 0.5 as we have chosen in our calcula- 
tions (all energy scales are normalized by ojafm in this paper) 
in this particular case. This is a rather artificial result of our 
model; nevertheless, it is an interesting behavior of the (D+S) 
type gap equations of Eq.(7-8). 

The summary for the (D+S) type gap is: (1) The mixed 
type gap solution is possible if the phonon interaction has a 
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FIG. 2: (Color online) (D+S) case. Calculated magnitudes of d-wave 
components (Ad/wafm) and s-wave components (As/u)afm) 
for different strength of phonon interaction N(0)Vp = 0.7,0.9, 
and 1.1. For all cases N(0)V M = 2.0, </>an = n/4 and 

Ulph/UAFM=0.5. 



FIG. 3: (Color online) (D+iS) case. Calculated magnitudes 
of d-wave components (Ad /loafm) and s-wave components 
(As/uiafm) for different anisotropy angles <f)AN = 7r/4,7r/2, 
and 7T. For all cases N(0)V M = 2.0, AT(0)Vp = 0.5 and 

UJph/u)AFM =0.5. 



proper anisotropy and its coupling strength is subdominant to 
the d-wave pairing strength; (2) T c is not enhanced as far as 
the d-wave gap remains finite regardless of the presence and 
magnitude of the s-wave gap component. 

B. (D+iS) case 

In this case, we assume AAFAi(k) — AdCOs(2</>) and 
A p h = iA a . Accordingly, A t (4>) — AdCOs(2</>) + iA s and 
E(<j>) = -\/£ 2 + A 2 d (4>) + A~J. The gap equations are 

<f>' 

iA a f» = -J2V P h(<t>-<f>')M<t>')x(<p',Uph) (12) 
<t>' 

The main difference from the (D+S) case is that because 
of the phase i between Ad(<fr) and A s , two gaps are only 
indirectly coupled through the quasiparticle energy E(<j>) in 
the pair susceptibility xi&i w A/,p/i)- Therefore the coupling of 
two gaps is much smoother than the (D+S) case. One conse- 
quence of it is that the gap equations can allow two transition 
temperatures T Cj d and T c ,s for each gap Ad and As, respec- 
tively. This feature is shown in Fig. 3 where N(0)Vm = 2.0 
(Xafm,d = 0.332) and N(0)V P = 0.55 are fixed and the 
anisotropy angle of phonon interaction 4>an is varied as 7r/4, 
7r/2, and tt (correspondingly \ ph)S = 0.243, 0.414 and 0.55). 



Main features are summarized as: (1) when X p h,s is too 
weak, the effect of phonon interaction is totally ignored and 
the gap and T c is solely determined by the AFM interaction 
(see 4>an = 7!"/4 case). (2) when X p h,s is much stronger 
than Xafm.d, the d-wave gap is totally suppressed and the 
gap and T c is determined solely by phonon interaction (see 
(f>AN — 7T case). This behavior is similar to the (D+S) case. 
(3) when the effective coupling strengthes are comparable, 
namely \afm,d ~ X p h,s< both d-wave and s-wave gaps have 
separate transition temperatures T c £> and T c g and two gaps 
can coexist at low temperatures (see 4>an — k/2 case). In 
this case, T Ci r> is always higher transition temperature and this 
transition temperature is just the same as the pure d-wave T c o 
without phonon interaction. At the lower transition tempera- 
ture T c ,s, the s-wave gap iAs appears and the magnitude of 
d-wave gap An accordingly decreases. 

To clarify the roles of the anisotropy and the strength 
of phonon interaction, in Fig.4 we fixed anisotropy an- 
gle as (f>AN = t> which simulates a perfectly isotropic 
phonon interaction, and varied the interaction strength. Fig.4 
shows the results of these calculations. With N(0)Vm =2.0 
(A AFM =0.332), we varied N(0)V P (= 0.3, 0.4, 0.5; accord- 
ingly A s =0.3, 0.4, 0.5). This plot demonstrates that the 
anisotropy of phonon interaction doesn't play a particular role 
in the case of the (D+iS) type gap. The mixed gap solution 
of the (D+ iS) type is still possible when Xafm,d and X p h,s 
are of the comparable strength regardless of the anisotropy of 
phonon interaction. 
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FIG. 4: (Color online) (D+iS) case. Calculated magnitudes of d- 
wave component (Ad/iOAFKi) and s-wave component (A s /uafm) 
for different strength of phonon interaction N(0)Vp — 0.3, 0.4, and 
0.5. For all cases N(0)Vm = 2.0, <f>AN = ir, uj ph /ujAFM=0.5. 



FIG. 5: (Color online) (DxFAf+Dph) case. Calculated mag- 
nitudes of T>afm component (Adi/uAFNi) and D p h compo- 
nent (A^/^afm) for different anisotropy angles 4>an = 
7r/8,7r/4,7r/2,and7r. For all cases N(0)V M = 2.0, N(0)V P = 0.5 
and Lo ph /ijJAFM=0.5. 



The summary for the (D+iS) type gap is: (1) The mixed 
type gap solution is possible if the phonon interaction \ p h,s 
is comparable but still subdominant to the AFM interaction 
strength Xafm,d- (2) In contrast to the (D+S) case, separate 
two transition temperatures T c .d and T c .s exist for each gap 
Art and As; T c ,d is always higher than T c .s and only below 
T c 5 two gaps coexist. (3) As in the (D+S) case, T c is not 
enhanced as far as the d-wave gap remains finite regardless 
of the presence and magnitude of the s-wave gap component. 
(4) Anisotropy of the phonon interaction is not particularly 
necessary. 



C. (D A FM+Dph) case 

This is the case that phonon interaction also supports 
the d-wave pairing. Because it is clear that a perfectly 
isotropic phonon interaction has null effect on the d-wave gap, 
anisotropy of the phonon interaction is crucial for this case. 
The total gap is, therefore, A t ((f>) = (A^i + Ad2)cos(20) 
and the gap equations are written as 
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The above two gap equations can be combined to a single 



FIG. 6: (Color online) (DAFM+D P h) case. Calculated magnitudes 
of the total d-wave gap (Ad-tot /^afm) for different strength of 
phonon interaction N(0)Vp = 0.0,0.2,0.4, and 0.6. For all cases 
JV(0)Vm = 2.0, 4>an = tt/4 andui ph /ujAFM=0.5. 
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gap equation but we keep the separate form in order to trace 
the effect of the phonon interaction. In contrast to the previous 
cases, T c can be dramatically enhanced when the phonon in- 
teraction Vph(A<t>) has a proper degree of anisotropy. Other 
authors 10 also obtained a similar result using different ap- 
proaches. In Fig. 5, we plot A^i and Ad2 separately; between 
the same symbols the smaller value is the phonon induced d- 
wave gap and the larger one is the AFM induced d-wave 
gap Adi- We change the anisotropy angle of V p h(A(j)) with 
the fixed interaction strengths of N(0)Vm = 2.0 (Xafm,d — 
0.332) and N(0)V p — 0.5. For the anisotropy angle 4>an — 
7r/2 and ir, the phonon interaction has absolutely no effect 
on the d-wave pairing; this is simply because the d-wave 
projected average Eq.(lO) becomes zero for these two com- 
mensurate angles with our model potential Eq.(3). With a 
stronger anisotropy (c/)an — tt/4; \ P h.D = 0.12), the phonon 
scattering sees only a limited part of the d-wave gap, there- 
fore the d-wave gap can be seen effectively as a s-wave gap 
for the anisotropic phonon. For a even stronger anisotropy 
(4>an — n/8;\ p h,D — 0.098), too narrow scattering angle 
reduces the effective coupling strength. In our model po- 
tential, 4>AN = tt/4 is the optimal anisotropy. In fact, even 
weaker anisotropy of 7r/4 < 4>a n < vr - except the special 
commensurate angles (f>AN — tt/2 and ir - produces some 
uncompensated effective interaction \ p h.D and boosts the d- 
wave pairing. This is an artifact of the model potential Eq.(lO) 
which always emphasizes the small angle scattering. There- 
fore, for some real phonons, it is possible to have a repulsive 
(destructive for d-wave pairing) X p h.D if large angle (around 
Acj) — tt/2) scattering overweighs small angle scattering. 

Fig. 6 shows the results of total gap Ad-tot — Adi + Ad2 
with the optimal anisotropy angle 4>an — 7r /4 and vary- 
ing phonon coupling strength iV(0)Vp=0.0, 0.2, 0.4, and 0.6 
(\ p , hD = 0.0, 0.047, 0.095, and 0.142, respectively). The re- 
sults demonstrate that the relatively weak phonon interaction 
Xph.D can significantly increase T c from the pure AFM inter- 
action induced T c o. 

Considering the experimental findings about phonons in 
the high-T c cuprates, the best candidate for the anisotropic 
pairing phonon is Bi 9 buckling mode of the plane oxygen 
motion 2 . Devereaux and coworkers extensively analyzed 
the behaviors of Bi g mode and found : (1) Bi ff mode has 
a strong anisotropic coupling matrix element concentrating 
around antinodal points; (2) the maximum coupling strength 
can reach as high as A ~ 3. Compared to our model calcu- 
lations, the optimal anisotropy for the d-wave pairing 4>an — 
7r/4 is just about the same degree of anisotropy of the Bi g 
mode. Regarding the coupling strength, for example, in Fig. 6 
the maximum coupling strength X max is N(0)V p =0.6 and the 
corresponding average strength is \ p h,D= 0.142. These val- 
ues are much smaller than the ones extracted by Devereaux 
and coworkers 2 ^. Larger values of coupling constants would 
be more realistic. However since our study is limited within 
the BCS gap equations we chose smaller values of coupling 
constants. With larger coupling strength (A > 0(1)), we 
should go for the strong coupling theory for the reliable anal- 
ysis. Also for the same phonon, the coupling strength ex- 
tracted from ARPES experiments and the coupling strength 
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FIG. 7: (Color online) (Pafm+Oph) gap- T c (symbols) calculated 
as a function of lo v h /luafm for different strength of phonon interac- 
tion N(0)Vp = 0.0, 0.2, 0.4, and 0.6, normalized by T c0 the transi- 
tion temperature with N(0)Vp = 0.0. For all cases, N(0)V M = 2.0 
and 4>an = it/ 4. Lines are the results of the analytic formula of T c 
Eq.(15). The symbols and lines of the same colors have the same 
parameters set. 



for the pairing problem can be very different due to the on- 
site Coulomb interaction. 

Another important anisotropic phonon in Bi-2212 is the 
breathing mode which strongly interacts near the nodes. The 
effect of this phonon is partially included in our calculations 
because our model phonon interaction Eq.(3) only constraints 
the scattering angle by 4>an but doesn't constraint the incom- 
ing k and outgoing momenta k (</> and cj) in our model). The 
anisotropic phonon scattering concentrating near nodes will 
be mostly averaged out by the sign changing d-wave gap and 
will be of no importance for the d-wave pairing. For realistic 
estimate of the phonon effects in HTC compounds, all impor- 
tant phonons should be included to calculate X p h.d using full 
coupling matrix elements. This is beyond the scope of the 
current paper. 

Since we found that an anisotropic phonon can dramat- 
ically enhance T c of d-wave pairing, we have to consider 
the isotope effect of the phonon, which is reported to be 
anomalously small, in particular, near optimal doping re- 
gion by many experiments' 2 . For this purpose, Fig. 7 shows 
the calculated T c 's (symbols) from Eq. (13-14) as a function 
of the phonon energy cutoff 0J p h/^AFM for various phonon 
coupling strength N(0)Vp = 0.0,0.2,0.4 and 0.6, respec- 
tively. The BCS theory predicts T c ~ u p h (a = 0.5), 
which should show up as a linear lines in Fig. 7. Our nu- 
merical results show much weaker power than the linear one 
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indicating a < 0.5. For more analytic investigation for 
the isotope effect and the origin of the T c enhancement by 
phonon in the (PAFM+^ph) case, we derived an analytic 
T c equation. The gap equations Eq.(13-14) can reduce to 
a single T c equation by adding two equations and taking a 
limit of A d i,A d2 — > 0. The pair susceptibility xi^AFM^h) 
(Eq.(9)) is integrated out using BCS approximation which is 
fv D d£[tanh « ln[2C lu d /ttT] (C = e< w 1.7807), 

valid when ojo/2T c 3> O(l). Finally we obtained the T c 
formula of the (Daf M+^pk) case as 



T c ~ 1.13 w 



AFM 



u ph 



-1/A t 



(15) 



where A t = (Xafm + X ph ), Xafm = Xafm/Xi and 
X p h = Aph/At, and Xafm and \ v h are the dimensionless ef- 
fective coupling constants obtained from Eq.(lO) with d-wave 
gap average for both couplings, loafm and uj p h are the en- 
ergy cutoffs of the AFM interaction and phonon interaction, 
respectively. As mentioned above, this T c formula is derived 
assuming ujM,ph/2T c ^> 0(1) and this condition is well satis- 
fied when At < 0(1). In Fig. 7, we also plot the results of T c 's 
from Eq.(15) (lines) with the same parameters and compare 
them with the numerical calculations (symbols) from Eq.(13- 
14). It shows a reasonably good agreement between two re- 
sults although the deviations increase with increasing the cou- 
pling constant N(0)V p (or \ p h,D) as expected. However, the 
overall and quantitative behavior of T c is well captured by 
Eq.(15). 

Now we are in position to read the phonon isotope coeffi- 
cient a from Eq.(15), which is 



1 



A 



ph 



2 ph 2 Xafm + X p h 



(16) 



For example, with representative values of Xafm — 0.33, 
Xph = 0.1 and uj p h/uJAFM — 0.5, we obtain a 0.116, 
which is pretty small value compared to the standard BCS 
value of 0.5 while T c is enhanced by about 100% from T c o the 
transition temperature without phonon interaction (see Fig. 7; 
N(0)V P = 0.6 is A p/l = 0.142). This isotope coefficient 
equation Eq.(16) provides a very plausible resolution why the 
isotope coefficient is so low near optimal doping region where 
T c is the highest and increases toward underdoped regime (de- 
creasing T c ). The phonon coupling strength (X p h) is likely un- 
changed with doping but the effective coupling strength of the 
AFM mediated interaction (Xafm) will change sensitively 
with doping. If we assume that Xafm increases with increas- 
ing doping but X p h remains constant, Eq.(16) describes the 
general trend of the experimentally observed oxygen isotope 
effect acPl 



D. (DAFM+Dph +iS ph ) case 
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FIG. 8: (Color online) (Dafm+D p /i +iS P h) case. Calculated mag- 
nitudes of total d-wave component (Ad = Aafm + A p h) and 
s-wave component (A s ) for different strength of phonon interac- 
tion N(0)Vp2 = 1.0,1.1 and 1.2. For all cases N(0)V M = 2.0, 
iV(0)Vpi = 0.6, 4>an = 7r/4, and lo p /luafm=0.5. 



(13), and (14) and the total gap will be A t (<j)) = (A dl + 
A d2 ) cos(2</>) + iA s . In this case, phonon interaction(s) me- 
diate both d-wave and s-wave pairings, which seems contra- 
dicting to a common knowledge. The fact is that phonon alone 
- no matter how anisotropic - cannot induce a d-wave gap but 
it can boost it if a d-wave gap is formed by other interaction. 
Then as we found in the (D+ iS) section, a phonon interaction 
can add iS component at lower temperatures. The interesting 
point is that the phonon interaction boosting d-wave gap and 
the phonon interaction inducing iS gap can originate from the 
same phonon if the anisotropy and the interaction strength are 
properly tuned. Of course more general case is that two differ- 
ent phonon modes play separate roles, respectively. For clear- 
ness, we rewrite the gap equations describing (Dafm+^ p Ii +i 
S p h) case below. 

A dl ((j)) = -J2Vafm{<!>- <l>')M<t>'M<l>' ,uafmUD 

A d2 (^) = -J2vi, P h(<fi-<i>')M<i>'M<f>\upi), (is) 
<p' 

iA s (<t>) - -5^2,p/,(0-^)A t r>')x(0\wp2) (19) 



Combining the {T>AFM+^ph) and (D+ iS) case, we can 
consider the (Dafm+^pH +*S p / l ) type solution. This is in- 
deed a natural solution of the combined gap equations (12), 



where the pair susceptibilities xi&i are defined in Eq.(9) 
and the quasiparticle energy is E((j)) — -y/£ 2 + A^(0) + 
with A d {<t>) — (Adi + A d 2) cos(20). For generality we 
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FIG. 9: (Color online) (Dafm+D p j, +i S p /i) case. Normalized spe- 
cific heats C(T)/C(T C ) for different strength of phonon interaction 
N(0)Vp2- There are two transition temperatures T c i (higher) and 
T C 2 (lower). Same parameters as in Fig. 8. 



assume two different phonon interactions Vi, p h and V^p/t, 
respectively but they can be the same phonon as explained 
above. 

In Fig. 8 we plot the solutions of A^ = A^i + Ad2 
and A s components separately. We fix the AFM interac- 
tion N(0)Vm=2.Q and the d-wave boosting phonon interac- 
tion N(0)Vpi=0.6 as a case of Fig. 6 and vary the phonon 
interaction iV(0)Vp2 inducing iS gap component as 1.0, 1.1, 
and 1.2, respectively. T c is enhanced by N(0)Vpi and doesn't 
change for different N(0) Vp% values because it is determined 
only by the d-wave gap component. We also chose the same 
anisotropy (f>AN = 7r /4 for both phonons Vi >p h and Vi,pK f° r 
convenience. Anisotropy condition for the d-wave boosting 
phonon is important; optimal anisotropy is 4>an — tt/4 for 
our model. However, the phonon interaction Vi, v h inducing 
iS gap doesn't need to be anisotropic. 

This type of gap solution can provide a possible resolution 
to the recent experiments of HTC which indicate the presence 
of a s-wave component at low temperatures. The anomalous 
tunnelling conductance in YBCC 3 , non-mono tonic behavior 
of the penetration depth at low temperatures and its isotope 
effect 4 , etc indicate a mixed gap of the (D+iS) type at low 
temperatures and a phonon effect with it. Fig. 9 show the nor- 
malized specific heat C(T)/C(T c i) for the gap solutions of 
Fig. 8. General feature is that it shows a d-wave type specific 
heat below T c i and then turns into a s-wave type behavior 
below T C 2 where iS gap component opens. Since the sec- 
ond transition at T C 2 is a true second order phase transition, 
the specific heat exhibits a jump. But the size of jump is far 
smaller than the BCS value of AC = 1.43 C(T c2 ). The rea- 
son is because when iS gap opens at T C 2, the magnitude of 
d-wave gap also abruptly decreases so that the amount of the 
specific heat jump is partially cancelled. Nevertheless, careful 
experiments should be able to discern the jump in the specific 



FIG. 10: (Color online) (Dafm+O v k +i S p h) case. Normalized su- 
perfluidity density Ai(0)/A|_(T) for different strength of phonon in- 
teraction N(0)Vp2- Same parameters as in Fig. 8. 



heat at the second transition temperature T C 2 if this scenario 
is true. 

Fig. 10 shows the normalized superfluid density 
A|(0)/A|(T) for the gap solutions of Fig.8. The gen- 
eral features are the same as in the specific heat: d-wave 
behavior in between T cl and T c2 , and s-wave behavior 
below T C 2. In particular, the additional sharp increase of 
XI 2 (T) below T c2 has some resemblance to the experimental 
observations 4 although the flattened out behavior at very low 
temperatures is not clearly observed experimentally. 

We propose a firm experimental test for this scenario. The 
d-wave boosting phonon interaction does not have a minimum 
interaction strength; as far as it has a proper anisotropy it will 
enhance the d-wave pairing as much as with its strength. For 
the iS component to appear below T c2 the phonon interaction 
needs a minimum strength to overcome a part of d-wave com- 
ponent. Therefore, a controlled suppression of d-wave gap 
(with non-magnetic impurities not to weaken the s-wave pair- 
ing) will widen a window of low temperature region where the 
i S component pops up. 



m. CONCLUSIONS 

In this paper we studied the effects of phonon interaction 
on the superconducting pairing in the background of d-wave 
gap already formed by the AFM interaction. In particular, we 
studied the role of anisotropy of the phonon interaction and 
possible multigap type solutions within a generalized BCS 
theory. For many cases the anisotropy of the phonon inter- 
action is a crucial condition for interesting interplay with the 
AFM interaction but not always; the (D+S) type gap need 
anisotropic phonon but (D+iS) type gap is more tolerable with 
the anisotropy of phonon interaction because two gaps are 
more gently coupled in the (D+iS) case. In both cases, T c is 
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not enhanced at all by phonon interaction as far as the d-wave 
gap component remains finite. 

Anisotropic phonon can boost T c together with the AFM in- 
teraction in the {Dafm+D p }i) type solution. With numerical 
calculations and analytic T c equation, we showed that T c is 
enhanced dramatically following the modified BCS exponen- 
tial form T c ~ exp ( — 1/A t ) with the total coupling strength 
At = Xafm + X p h- This type of solution can also explain 
how and why the phonon isotope effect is strongly reduced 
despite the large enhancement of T c by phonon. Our isotope 
coefficient formula Eq.(16) provides a good description of the 
overall trend of oxygen isotope coefficient acr^. 

Also a combined type solution of (Dafm + D p h + iS) gap 
is considered. This type of gap not only shows the features 
of the (Dafm + -Dph) g a P - enhanced T c and small isotope 
effect - but also shows the appearance of a s-wave component 
at low temperatures. The latter feature can provide a micro- 
scopic origin of the small s-wave component suggested from 
recent experiments of tunnelling conductance 3 and penetra- 
tion depth measurements 4 . Calculations of the superfluidity 
density showed a qualitative agreement with experiments of 
the abrupt increase of A^ 2 (T) at lower temperatures 4 . How- 
ever, this scenario should be confirmed by experimental ob- 
servation of a specific heat jump at the second low transition 



temperature T C 2 for the s-wave gap component as shown in 
Fig.9. We proposed a systematic impurity doping study to test 
this scenario. 

Finally, we would like to restate the limitation of our work. 
First, the dynamical interplay between phonon(s) and AFM 
interaction is not properly treated. This problem was already 
studied by many authors^2LL!J We bypassed this complicated 
problem on the phenomenological basis; after all correlation 
effects taken place we assumed that fermion quasiparticles re- 
main, and AFM spin fluctuations and several phonon modes 
are considered as experimentally defined quantities after mu- 
tual screening and renormalization. Second, when we con- 
sider the s-wave pairing, we should have included the on-site 
Coulomb interaction in the HTC cuprates in addition to the 
phonon interaction. Therefore, X p h,s m this paper should be 
considered as a renormalized quantity after including the on- 
site Coulomb interaction. Perhaps, this is the reason why our 
Xph s (~ 0-1) is order of magnitude smaller than the experi- 
mental values^ \ exp ~ 0(1). 

I thank J. Zaanen for drawing my attention to this problem 
and discussions during MSM07 in Khiva. Y. B. was supported 
by the KOSEF through the CSCMR and the Grant No. KRF- 
2007-070-C00044, KRF-2007-52 1 -C0008 1 . 
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